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We review the 0(a) improvement of lattice QCD with special emphasis on the motivation for performing 
the improvement programme non-perturbatively and the general concepts of on-shell improvement. The present 
status of the calculations of various improvement coefficients (perturbative and non-perturbative) is reviewed, as 
well as the computation of the isospin current normalization constants Zk and Z\. We comment on recent results 
for hadronic observables obtained in the improved theory. 



1. INTRODUCTION 

The leading cutoff effects in lattice QCD with 
Wilson quarks are proportional to the lattice 
spacing a and can be rather large for typical val- 
ues of the Monte Carlo (MC) simulation param- 
eters . Decreasing the lattice spacing at con- 
stant physical length scales means larger lattices 
and therefore rapidly increasing costs of the MC 
simulations. 

An alternative method to reduce cutoff effects 
in lattice field theories is due to Symanzik 
He has shown that the approach of Green's func- 
tions to their continuum limit can be accelerated 
by using an improved action and improved (com- 
posite) fields. A considerable simplification is 
achieved if the improved continuum approach is 
only required for on-shell quantities such as parti- 
cle masses and matrix elements of improved fields 
between physical states j6|j^. 

Symanzik improvement may be viewed as an 
extension of the renormalization programme to 
the level of irrelevant operators. While the 
structure of the possible counterterms is dic- 
tated by the symmetries, their coefficients have 
to be fixed by appropriate improvement con- 
ditions. Although they can be estimated in 
(tadpole-improved |^) perturbation theory, a 
non-perturbative determination of the improve- 
ment coefficients through MC simulations is 
clearly preferable. 

To 0(a), we have recently carried out the 
non-perturbative improvement programme in 

* Invited talk presented at the workshop Lattice QCD on 
parallel Computers, CCP, Tsukuba, Japan, March 1997. 



quenched lattice QCD, thus leaving residual cut- 
off effects of O(a^) only. The details of the cal- 
culations are given elsewhere In this re- 
port we emphasize the general concepts and give 
a short account of the main results. We also dis- 
cuss progress that is currently made in applying 
the improvement programme to full QCD with 
A^f = 2 flavours 

We start by outlining the importance of re- 
moving the linear (in a) lattice artefacts non- 
perturbatively (sect. 2). We then review on-shell 
improvement in the framework of Symanzik's lo- 
cal effective theory (sect. 3) and the general con- 
cept of non-perturbative improvement (sect. 4). 
The practical calculations are done in finite space- 
time volume with Schrodinger functional bound- 
ary conditions which we introduce in sect. 5. We 
then review the computations of the various im- 
provement coefficients and the current normaliza- 
tion constants in quenched lattice QCD. Resid- 
ual lattice artefacts are discussed in sect. 9 and 
the status of 0(a) improvement in full QCD in 
sect. 10. 

2. WHY NON-PERTURBATIVE 0(a) IM- 
PROVEMENT? 

In this workshop two main alternative ap- 
proaches to reducing the cutoff effects of lattice 
QCD are discussed [ p"5|Jl6| . Before giving some 
detailed motivation for what is discussed here, 
let us point out the relation to these talks. The 
goal of the perfect action approach [|l5| is the im- 
provement of the lattice theory to all orders of the 
lattice spacing. In determining the perfect action. 
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one - in practice - performs an expansion in pow- 
ers of the coupling constant. In contrast, the non- 
perturbative Symanzik improvement as discussed 
here proceeds non-perturbatively in the coupUng 
and order by order in the lattice spacing. On 
the other hand, "mean field improvement" fl^ 
is Symanzik improvement with improvement co- 
efficients approximated by mean field-improved 
perturbation theory l^. 

Let us explain in more detail the motivation for 
non-perturbatively removing the lattice artefacts 
that are linear in a. The first motivation is that 
these are large numerically. We give two examples 
to illustrate this point. 

The first example is a perturbative one. Con- 
sider the renormalized coupling g{L) defined in 
the Schrodinger functional scheme | p7| , plp9| ] for 
massless fermions. Its evolution from length scale 
L to scale 2L defines the step scaling function 
E The one-loop contribution of A^f Wil- 

son fermions to E has been calculated in ref. |l9| . 
We denote it by g^TVfSi^i . Its continuum hmit 
is given by the fermion contribution to the lowest 
term in the /3-function: 



a/L^o g Ni 



(1) 
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Without O(o)-improvement, Ei_i depends very 
strongly on a/L and shows a smooth approach 
to the continuum only for very large lattices (cf. 
fig. |). If one had MC data with L/a < 10 and 
such an error term, a naive continuum extrapola- 
tion of such data would overestimate the fermion 
contribution by 0(60%). On the other hand, af- 
ter 0(a)-improvement the lattice artefacts in this 
quantity are tiny. 

The second example is the current quark mass 
m defined by the PCAC relation. As we will dis- 
cuss in detail below, its value is independent of 
kinematical variables such as the boundary con- 
ditions. Dependences on such variables are pure 
lattice artefacts. We examined the current quark 
mass in the valence approximation by numerical 
Monte Carlo simulations and found |^ large lat- 
tice artefacts even for quite small lattice spacings 
(cf. fig. I). 
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Figure 1. One-loop contribution of a massless 
Wilson quark to the step scaling function S nor- 
malized to its continuum value. Open circles are 
for the Wilson action. The dotted curve describes 
the leading lattice artefact for small a/L, which 
is of the form log{a/ L)a/ L, while the dashed 
curve contains the pure a/i-term in addition. 
Filled circles illustrate the behaviour after 0(a)- 
improvement. 



Since lattice artefacts can be large numerically, 
it is desirable to remove them non-perturbatively 
in order not to leave over a^a or truly non- 
perturbative terms that may be noticeable. 

A second motivation for performing improve- 
ment non-perturbatively is closely related. Af- 
ter improvement one wants to compute physical 
quantities P, such as ratios of hadron masses for 
a range of lattice spacings and extrapolate them 
to the continuum limit. 

For example, assume that MC results P{a) are 
available for four points = Oq, jO^, ^Oq, jO^ 
and with statistical errors AP — A that are in- 
dependent of a. An extrapolation by means of a 
fit P{a) = Pc + 320? will yield an estimate for the 
continuum limit of P with error APc = 1.22 A. 
However, if the calculation was performed with 
approximate values of the improvement coeffi- 
cients, a small (?) linear error term 6 x a will be 
left in the results P{a). If ignored, it causes a sys- 
tematic error of size 0.35 6 oq in Pc- As an alter- 
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Figure 2. Dependence of curcnt quark mass m 
on the boundary condition and the time coordi- 
nate. The calculation is done on a 16 x 8 lattice 
a.t /3 — 6.4, which corresponds to a lattice spac- 
ing of a « 0.05 fm. "Boundary values" refer to 
the gauge field boundary conditions in the SF ; 
see sects. for details. 



native, one may, of course, perform a fit with the 
ansatz P{a) = Pc + Sia -I- S2a^. This gives a con- 
tinuum value with statistical error AP^ = 9.7 A. 

Trivial as it is, the above example illus- 
trates the important point: in order to perform 
a reliable continuum extrapolation, one must 
have theoretical information on the expected a- 
dependence. This is contained in the structure 
of Symanzik's effective action, which justifies a 
power series in the lattice spacing for the cor- 
rection terms. (For the purpose of extrapolating 
MC-data, the logarithmic dependence on a can be 
ignored). Under the same assumption, namely 
that lattice artefacts higher than can be ne- 
glected, the continuum extrapolation in the 0(a)- 
improved theory gives errors that are smaller by a 
factor 8 than the extrapolation in the unimproved 
(or approximately improved) case. 

Natural as it may seem to continue the non- 
perturbative Symanzik improvement beyond the 
first order in the lattice spacing, this appears im- 
possible in practice. The reason is that at or- 
der O(a^) a large number of improvement terms 



would have to be determined by MC calculations. 

3. ON-SHELL IMPROVEMENT 

3.1. Lattice QCD with Wilson quarks 

In this section we consider QCD on an infinitely 
extended lattice with two degenerate light Wilson 
quarks of bare mass mo ||]. The action is the 
sum of the usual Wilson plaquette action and the 
quark action 

5F[f/,?,V'] =a'5]^(^)(^ + "^o)^(^)' (2) 

X 

where a denotes the lattice spacing. The Wilson- 
Dirac operator 

3 

^=5E[(^^ + ^m)7m-«V*V^], (3) 

contains the lattice covariant forward and back- 
ward derivatives, and V*. The last term in 
eq. (^ eliminates the unwanted doubler states but 
also breaks chiral symmetry. As a consequence, 
both additive and multiplicative renornialization 
of the quark mass are necessary, i.e. any renor- 
malized quark mass mR is of the form 

mn = Zi„mq, Wq = toq - mc (4) 

where rric is the so-called critical quark mass. 

Chiral symmetry violation is more directly seen 
by studying the conservation of the isovector axial 
current A°;. The current and the associated axial 
density on the lattice are defined through 



A^(x) = ^(a:)7^75— V;(x), (5) 

P"(x) = ^(a:)j,^i>{x), (6) 

where r" are Pauli matrices acting on the flavour 
index of the quark field. The PCAC relation 

9^ a; (a;) = 2mP«(a;) + 0(a), (7) 
d, = + (8) 



then includes an error term of order a, which can 
be rather large on the accessible lattices , as 
shown in fig. |[ Above, 9^ and 9jt denote for- 
ward and backward lattice derivatives and m is 
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the iinrenormalized current quark mass at scale 
1/a. 

The isospin symmetry remains unbroken on the 
lattice and there exists an associated conserved 
vector current. However, it is often advantageous 
to use the current which is strictly local, 



(9) 



The conservation of this current is then also vi- 
olated by cutoff effects, and a finite renormal- 
ization is required to ensure that the associated 
charge takes half-integral values. 

3.2. Symanzik's local effective theory 

Near the continuum limit the lattice theory can 
be described in terms of a local effective the- 
ory D, 



ScS ^ Sq + aSi +a S2 + 



(10) 



where Sq is the action of the continuum theory, 
defined e.g. on a lattice with spacing e a. The 
terms Sk, k — 1,2,..., are space-time integrals 
of Lagrangians Ck {x) . These are given as general 
linear combinations of local gauge-invariant com- 
posite fields which respect the exact symmetries 
of the lattice theory and have canonical dimen- 
sion A + k. We use the convention that explicit 
(non-negative) powers of the quark mass m are 
included in the dimension counting. A possible 
basis of fields for the lagrangian Ci (x) then reads 



O5 — m^'ip^p, 



(11) 



where F^i^ is the field tensor and cr^t^ = -^[jp,, 7i/]. 

When considering correlation functions of local 
gauge invariant fields the action is not the only 
source of cutoff effects. If 0(x) denotes such a lat- 
tice field (e.g. the axial density or the isospin cur- 
rents of subsect. ||.l), one expects the connected 
n-point function 

Gnixi, . . . , Xn) - {Z^r W^l) ■ ■ ■ '^(a;n))co„ (12) 



to have a well-defined continuum limit, pro- 
vided the renormalization constant is correctly 
tuned and the space-time arguments cci, . . . , a;„ 
are kept at a physical distance from each other. 

In the effective theory the renormalized lattice 
field Z,f,(f>(x) is represented by an effective field, 

(j>cs{x) = (poix) + a<j)i{x) + a'^(j>2{x) + (13) 

where the (t>k{x) are linear combinations of com- 
posite, local fields with the appropriate dimension 
and symmetries. For example, in the case of the 
axial current (|^), (pi is given as a linear combina- 
tion of the terms 

{Or)l = i^^T-j^iD^ + D^]^, (14) 

The convergence of Gn{xi, . . . , x„) to its contin- 
uum limit can now be studied in the effective the- 
ory, 

Gn{xi,...,Xn) = {(po{xi) . ■ ■ (pQ{Xn)) con 

-a I d^y {(po{xi) . . . (pa{xn)Ci{y))^^^ 



(15) 



{<Po{xi) ■ ■ ■ (plixk) ■ ■ ■ (po{Xn))c 



k=l 



where the expectation values on the right-hand 
side are to be taken in the continuum theory with 
action 5*0. 

3.3. Using the field equations 

For most applications, it is sufficient to com- 
pute on-shell quantities such as particle masses, 
S-matrix elements and correlation functions at 
space-time arguments, which are separated by a 
physical distance. It is then possible to make use 
of the field equations to reduce first the number 
of basis fields in the effective Lagrangian Ci and, 
in a second step, also in the 0(a) counterterm cpi 
of the effective composite fields. 

If one uses the field equations in the Lagrangian 
Ci under the space-time integral in eq. (15), the 
errors made are contact terms that arise when y 
comes close to one of the arguments xi, . . . , Xn- 
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Taking into account the dimensions and symme- 
tries, one easily verifies that these contact terms 
must have the same structure as the insertions of 
(pi in the last term of eq. (^5|). Using the field 
equations in £i therefore just means a redefini- 
tion of the coefficients in the counterterm . 

It turns out that one may eliminate two of the 
terms in cq. (pl|). A possible choice is to stay 
with the terms Oi, O3 and O5, which yields the 
effective continuum action for on-shell quantities 
to order a. Having made this choice one may 
apply the field equations once again to simplify 
the term (j>i in the effective field as well. In the 
example of the axial current it is then possible to 
eliminate the term Og in eq. ( p^ ) . 

3.4. Improved lattice action and fields 

The on-shell 0(a) improved lattice action is 
obtained by adding a counterterm to the unim- 
proved lattice action such that the action Si in 
the effective theory is cancelled in on-shell am- 
plitudes. This can be achieved by adding lattice 
representatives of the terms d , O3 and O5 to the 
unimproved lattice Lagrangian, with coefficients 
that are functions of the bare coupling go only. 
Leaving the discussion of suitable improvement 
conditions to sect. 4, we here note that the fields 
O3 and C5 already appear in the unimproved the- 
ory and thus merely lead to a reparametrization 
of the bare parameters go and mg. In the fol- 
lowing, we will not consider these terms any fur- 
ther. There relevance in connection with massless 
renormalization schemes is discussed in detail in 
ref. §. 

We choose the standard discretization Ff^i, of 
the field tensor ^ and add the improvement term 
to the Wilson-Dirac operator (||), 

m ^ 

With a properly chosen coefhcient Csw(ffo)i this 
yields the on-shell 0(a) improved lattice action 
which has first been proposed by Sheikholeslami 
and Wohlert 0. 

The perturbative expansion of Csw reads Csw = 
1 + c^J^gl + 0(g4), with 00 c}^ - 0.26590(7). 

The 0(a) improved isospin currents and the ax- 



ial density can be parametrized as follows, 
(Ar); - ZA(l + 6Aamq){A^ + acAa^F'^}, 
{V,,)l - Zv(l + 6vamq){^; + acyd^T^,], 
(Pr^ = ^p(l + fepamq)P^ (17) 
where 

and the other fields have been defined in sub- 
sect. 2.1. We have included the normalization 
constants Za,v,P5 which have to be fixed by ap- 
propriate normalization conditions (cf. sect.^. 
Again, the improvement coefficients 6a,v,p and 
CA,v are functions of go only. At tree level of per- 
turbation theory, they are given by 6a = &p = 
6v = 1 and CA = cv = Due to the 

efforts of Liischer, Sint and Weisz ]lO| , p^ these 
coefficients are now known to one-loop accuracy, 

CA = -0.005680(2) xg2cF + 0(g;^), 

cv = -0.01225(1) X .ggCp + 0(5^), 

6a = 1 + 0.11414(4) X5o2CF + 0(g^), 

6v = 1 + 0.11492(4) xg2CF + 0(g4), 

6p 1 + 0.11484(2) x.g2Cp + o(.g^). (18) 

Here Cp = 4/3. It is not known why 6a, 6v and 
6p are numerically so close to each other to this 
order of perturbation theory. 

Non-perturbative determinations of Csw, ca, cy 
and 6v will be discussed in sections 6-8. 

3.5. The PCAC relation 

We assume for the moment that on-shell 0(a) 
improvement has been fully implemented, i.e. the 
improvement coefficients are assigned their cor- 
rect values. If O denotes a renormalized on-shell 
0(a) improved field localized in a region not con- 
taining X, we thus expect that the PCAC relation 

9^((Ar)«(x) O) - 2mR((PR)'^(a;) O) (19) 

holds up to corrections of order . At this point 
we note that the field O need not be improved for 
this statement to be true. To see this we use again 
Symanzik's local effective theory and denote the 
0(a) correction term in Ocs by 0i. Eq. ( |l9| ) then 
receives an order a contribution 

aiidM^T^^ix) - 2mR(FR)'^(a:)} </)i), (20) 
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which is to be evaluated in the continuum theory. 
The PCAC relation holds exactly in this limit and 
the extra term ( po| ) thus vanishes. 

4. NON-PERTURBATIVE 
IMPROVEMENT 

We have explained above how the form of the 
improved action and the composite fields is de- 
termined by the symmetries of the lattice action. 
The coefficients of the different terms need to be 
fixed by suitable improvement conditions. One 
considers pure lattice artefacts, i.e. combinations 
of observables that are known to vanish in the 
continuum limit of the theory. Improvement con- 
ditions require these lattice artefacts to vanish, 
thus defining the values of the improvement coef- 
ficients as a function of the lattice spacing. 

In perturbation theory, lattice artefacts can 
be obtained from any (renormalized) quantity by 
subtracting its value in the continuum limit. The 
improvement coefficients are unique.^ 

Beyond perturbation theory, one wants to de- 
termine the improvement coefficients by MC cal- 
culations and it requires significant effort to take 
the continuum limit. It is therefore advantageous 
to use lattice artefacts that derive from a sym- 
metry of the continuum field theory that is not 
respected by the lattice regularization. One may 
require rotational invariance of the static poten- 
tial V{r), e.g. 

V{v = (2, 2, l)r/3) - V{r = (r, 0, 0)) = , 

or Lorentz invariance, 

[i?(p)]2-[i?(0)p-p2=0, 

for the momentum dependence of a one-particle 
energy E. 

For 0(a) improvement of QCD it is advanta- 
geous to use violations of the PCAC relation (|^), 
instead. PCAC can be used in the context of 
the Schrodinger functional (SF), where one has 

^One might think there is a scheme ambiguity as in the 
general renormahzation of a theory. The "scheme" has, 
however, already been chosen implicitly, when one writes 
down the original (unimproved) lattice theory. It changes, 
for instance, when a different discretization of the pure 
gauge action is chosen. 



a large flexibility to choose appropriate improve- 
ment conditions and can compute the improve- 
ment coefficients also for small values of the bare 
coupling goi making contact with their perturba- 
tive expansions. A further - and maybe the most 
significant - advantage of the SF in this context is 
the following. In the SF we may choose boundary 
conditions such that the classical solution, called 
the induced background field, has non- vanishing 
components [0. Remembering eq. (p^), we 
observe that correlation functions are then sensi- 
tive to the improvement coefficient Csw already at 
tree level of perturbation theory. In general this 
will be the case at higher orders only. This is the 
basis for a good sensitivity of the improvement 
conditions to Csw 

As a consequence of the freedom to choose im- 
provement conditions, the resulting values of im- 
provement coefficients such as Csw, ca depend on 
the exact choices made for the improvement con- 
ditions. The corresponding variation of Csw, ca 
is of order a. This variation changes the effects 
of order in physical observables computed af- 
ter improvement. In principle this is irrelevant at 
the level of 0(a) improvement. Nevertheless, one 
ought to choose improvement conditions where 
such terms have small coefficients. The improve- 
ment conditions derived from the SF can be stud- 
ied in perturbation theory. Such a study pro- 
vided essential criteria for our detailed choice of 
improvement conditions. 

In the following section we introduce the SF to 
prepare for the improvement conditions that are 
detailed in sects. 

5. THE SCHRODINGER FUNCTIONAL 

5.1. Definitions 

The space-time lattice is now taken to be a dis- 
cretized hyper-cylinder of length T and circum- 
ference L. In the spatial directions the quantum 
fields are L-periodic, whereas in the Euclidean 
time direction inhomogeneous Dirichlet bound- 
ary conditions are imposed as follows. The spa- 
tial components of the gauge field are required to 
satisfy 

U{x,k)l^^^ = eMaCk), Ck = i(f>/L, (21) 
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with (f) = diag(0i, 02, (^3), and an analogous 
boundary condition with C" is imposed at Xq — T. 

With the projectors P± = i(l±7o), the bound- 
ary conditions for the quark and antiquark fields 
read 

P+^U„=0 = P, P~^xo=T - p\ (22) 

V^P_U„=n = P, V^P+U„=T = p'. (23) 
The functional integral in this situation [[l7pif| , 



Z[C',p',p';C,p,p] 



(24) 



fields 



is known as the QCD Schrodinger func- 
tional (SF). Concerning the (unimproved) action 
S, we note that its gauge field part has the same 
form as in infinite volume. The quark action is 
again given by eq. (H) , provided one formally ex- 
tends the quark and antiquark fields to Euclidean 
times Xo < and a;o > T by "padding" with 
zeros However, we will use a slightly more 
general covariant lattice derivative, 

V^^{x) = - [\^,U{x, n)^/j{x + afi) - ^(x)] , (25) 
a'- ■' 

where Aq = 1 and = exp(ia0/_L). This modifi- 
cation of the covariant derivative is equivalent to 
demanding spatial periodicity of the quark fields 
up to the phase exp(z6'). We thus have the an- 
gle 9 as an additional parameter that plays a role 
in the improvement condition for the coefhcient 

CA0. 

We are now prepared to define the expectation 
values of any product O of fields by 



(O) 



1 



Oi 



fields 



(26) 



p'=p'=p=p=0 



Apart from the gauge field and the quark and 
anti-quark fields integrated over, O may involve 
the "boundary fields" at time xq = 0, 



C(x) 



C(x) 
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(27) 



and similarly the fields at xq — T. Note that 
the functional derivatives only act on the Boltz- 
mann factor, because the functional measure is 
independent of the boundary values of the fields. 



5.2. Continuum limit and improvement of 
the Schrodinger functional 

Based on the work of Symanzik (2^j2^ and 
explicit calculations to one-loop order of pertur- 
bation theory one expects that the SF 
is renormalized if the coupling constant and the 
quark masses are renormalized in the usual way 
and the quark boundary fields are scaled with 
a logarithmically divergent renormalization con- 
stant. 

As in the case of the infinite volume theory 
discussed in subsect.|[2, the cutoff dependence 
of the SF may be described by a local effective 
theory. An important difference is that the 0(a) 
effective action Si now includes a few terms local- 
ized at the space-time boundaries ||^ . Such terms 
then also appear in the 0(a) improved lattice ac- 
tion. However, by an argument similar to the 
one given at the end of subsect.||.5, it can be 
shown that they only contribute at order to 
the PCAC relation and the chiral Ward identity 
considered in sect. 5. In the calculations reported 
below, the inclusion of boundary counterterms is, 
therefore, not required. 

6. IMPROVEMENT OF ACTION AND 
AXIAL CURRENT 

We now proceed to sketch the non-perturbative 
calculation of the coefficient Cs„ 
Using the operator 



0'^ = a«5]C(y)75^C(z), 



we define the bare correlation functions 

/A(a;o) = -i(AS(.T)0'^), 
/p(.To) = -^,{P-ix)0'^). 



(28) 



(29) 



In terms of /a and /p the PCAC relation for the 
unrenormalized improved axial current and den- 
sity may be written in the form 



t^Q/A(a^o) + CAadodofpixo) 
2/p(xo) 



(30) 



We take this as the definition of the bare current 
quark mass m. The renormalized quark mass tor 
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appearing in eq. ( [l9| ) is then given by 



Tor 



m 



ZA{l + bAamq) ^Q^^2) 



Zp(l + bpamq) 

At fixed bare parameters, mpt and hence also the 
unrenormalized mass m should be independent 
of the kinematical parameters such as T, L and 
xq. This will be true up to corrections of order 

, provided Cgw and ca have been assigned their 
proper values. The coefficients may, therefore, be 
fixed by imposing the condition that m has ex- 
actly the same numerical value for three different 
choices of the kinematical parameters. 

For the rest of this section we set T — 2L, 9 = 

and {(j)i,(j}2,h) = B (-^'O''^) '('^i''^2''^3) = 
^ (— 57r, 2tt, Stt). An important practical criterion 
for choosing these particular boundary values has 
been that the induced background field should be 
weak on the scale of the lattice cutoff to avoid 
large lattice effects. On the other hand, the ef- 
fects of order a, which one intends to cancel by 
adjusting Cgw should not be too small: otherwise 
one would be unable to compute Csw accurately. 
The above boundary values represent a compro- 
mise where both criteria are fulfilled to a satisfac- 
tory degree on the accessible lattices. 

Another mass, m', may be defined in the same 
way by interchanging C and C . The PC AC rela- 
tion then implies that the mass difference m — m! 
is of order a? if the coefficients Csw and ca are ap- 
propriately chosen. Our intention in the following 
is to take this as a condition to fix Csw 

Before we are able to do so, we must eliminate 
the coefficient ca, which is not known either at 
this point. To this end first note that 

m(xo) = r(xo) + CAs(xo), (32) 
where r and s are defined through 

r{xo) = K^o + 5o)/a(xo)//p(xo), (33) 

(34) 



s(xo) = \adtdof-p{xQ)l fp{xQ) 

(for clarity the dependence on the time coordi- 
nates is now often indicated explicitly) . The other 
mass m' is similarly given in terms of two ratios 
r' and s'. It is then trivial to show that the com- 
bination 

M(xo,yo) ^ m[xo) - s(xa) — j — ^ 77 — r- (35) 
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' I I I I I I I I 

0.02 0.04 0.06 
aM 

Figure 3. Mass difference AM on a 16 x 8"^ lattice 
as a function of the quark mass M at /? = 6.4 and 
Csw = 1.777. 

aAM 




Figure 4. Determination of Csw at /3 = 6.4. The 
dashed line indicates the tree-level value AM^°^ 
appearing on the right-hand side of the improve- 
ment condition. 



is independent of ca, viz. 

M{xo,yo) = r{xQ)-s{xa)-— — — . (36) 

s{ya) - s'iyo) 

Furthermore, from eq. ( ^^ one infers that M 
coincides with m up to a small correction of or- 
der (in the improved theory); M may hence be 
taken as an alternative definition of an unrenor- 
malized current quark mass, the advantage being 
that we do not need to know ca to be able to 
calculate it. 

We now continue to discuss the condition that 
determines Csw • If we define M' in the same way 
as M, with the obvious replacements, it follows 
from the above that the difference 



AM = MRT,\T) 



M' 



4^ ' 4^ / 



(37) 
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must vanish, up to corrections of order o^, if Csw 
has the proper value. The coefhcient may hence 
be fixed by demanding AM = AM*^°^ Here 
AM^"^^ , the value of AM at tree-level of perturba- 
tion theory in the 0(a) improved theory, is cho- 
sen instead of zero, in order to cancel a small 
tree-level 0(a) effect in Csw This is more a mat- 
ter of aesthetics than of practical importance (cf. 
fig. ^). In order to complete the improvement con- 
dition, we further give the precise definition of the 
quark mass: we evaluate AM at M(iT, \T) = 0. 
Again this choice is inessential, changing Cgw only 
to O(aM), which are negligible effects indeed (cf. 
fig. I). 

At each of our nine values of go, we compute 
AM for M = for at least three values of Csw, 
and AM = AM*^°) is solved for Csw by a linear fit 
of AM as a function of Csw Representative data 
are displayed in fig.0. 

In the range < g^, < 1, the results for Csw are 
well represented by [|ll| 



1 - 0.656 - 0.152 g^ - 0.054 g^ 



(38) 



" 1 - 0.922 gl 

In fig. 1^ we compare our results to one- loop bare 
perturbation theory and also to tadpole-improved 
perturbation theory [p| , for which we have used 



l + (cW-l/4)5i 



P ! 



(39) 



where gl = g'^/i 



Here Uq is the average 



plaquette in infinite volume. 



CA 



In a similar way In) we obtained (0 < g§ < 1) 
1 - 0.748 g2 



-0.00756 X g^- 



(40) 



1 - 0.977 5o'' 

as displayed in fig. ||. Both eq. (^ and eq. (40) 
deviate substantially from the one-loop result ex- 
cept for values of g^ as small as < 1/2. 

7. CURRENT NORMALIZATION 

7.1. Chiral Ward identities 

For zero quark masses, chiral symmetry is ex- 
pected to become exact in the continuum limit. It 
has therefore been proposed to fix the normaliza- 
tion constants Za and Zy of the isospin currents 
by imposing the continuum chiral Ward identities 
also at finite values of the cutoff CQl-lESl. 




Figure 5. Non-perturbative improvement coef- 
ficient Csw. Results from one- loop bare and tad- 
pole improved perturbation theory are denoted 
by dotted lines and crosses, respectively. 



In the case of the axial current the relevant 
Ward identity can be written in the form 

da, {x)e^'^ {A; {x)AUy)Q^)^ 2r{V^ {y)Q^), (41) 

dR 

where the integral is taken over the boundary of 
the space-time region R containing the point y 
and Q"" is a source located outside R. In view 
of on-shcU improvement, it is important to note 
that all space-time arguments in eq. (41) are at 
non-zero distance from one another. 

For the source field in eq. (^Tj) we choose Q° = 

defined similarly using the primed fields. The re- 
gion R is taken to be i? = {x, ti < xq < ^2}- 

We define the following correlation functions, 
using the unrenormalized improved currents 
(Ai)- ^ A-+ acAd^P- and (Vi)" - + 



/aa (3^0,2/0) 



bc^cde 
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Figure 6. Non-perturbative improvement coeffi- 
cient ca- Results from one- loop bare and tadpole 
improved perturbation theory are denoted by dot- 
ted lines and crosses, respectively. 



x(0"^(Ai)g(x)(Ai)|5(y)0'=), (42) 



^^{0'%V{)l{x)0^), (43) 



(44) 



At TTT-q = 0, and using the correct values of 
Csw and Ca, a lattice version of the chiral Ward 
identity (|l|) is 

^i/iA(a;o,yo)-/i + 0(a2), 



xo> Va- (45) 



Compared to eq. ( [4l|) we have set v — Q and in- 
cluded an additional summation over y, thus ob- 
taining the isospin charge. In deriving eq. (^) we 
have used the fact that the action of the latter 
on the chosen matrix elements can be evaluated 
thanks to the exact isospin symmetry on the lat- 
tice. We further exploited the conservation of the 
axial charge (at zero mass). 

Since the isospin symmetry remains unbroken 
for non-zero quark mass, one need not restrict 
the normalization condition of the vector current 



to the case mq = 0. Using arguments similar to 
those in the derivation of eq. (|45|), one obtains 



Zv(l + &vamq)/^(a;o) = /i + 0(a2). 



(46) 



Note that the improvement coefficient cy is not 
needed here, because the tensor density does not 
contribute to the isospin charge. 

7.2. Normalization conditions 

Starting from eqs. (^,|4^) we impose the fol- 
lowing normalization conditions on the axial and 
vector currents at vanishing quark mass, mq = 0, 

^A/AA(ir,ir) = A, T = 9L/4, (47) 
Zwfy{\T) = /i, T = 2L, (48) 

where we have set C = C = 9 = 0. In or- 
der to guarantee that cutoff effects of matrix el- 
ements of the renormalized currents vanish pro- 
portional to when approaching the continuum 
limit, we complete the normalization conditions 
eqs. ( ^ , ^8| ) by scaling L in units of a physical 
scale, 



L/a — 8 at go = 1, 
L/ro — constant. 



(49) 



Here rp ~ 0.5fm is derived from the force be- 
tween static quarks as explained in ref. [p2[. In 



the numerical simulations, the choices (t47-49) to 



gether with a specific definition of realized 
by smooth inter-/extrapolations 

Of course, the details of these choices are irrel- 
evant. Nevertheless, one must take care not to 
induce large O(a^) effects through an "unfortu- 
nate" choice. For this reason the Ward identities 
were studied in perturbation theory before fix- 
ing the above details. In addition we checked in 
the course of the simulations that Zy, Za do not 
change appreciably when parameters such as ( ^ ) 
or 9 are changed within reasonable limits. 

7.3. Results 

Our numerical results are again well repre- 
sented by rational functions 



Za 



1 - 0.7663 + 0.0488 

1 - 0.6369 g^ ' 

1 -0.84965^ + 0.06105^ 
1 - 0.7332 5^ ' 



(50) 
(51) 
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Figure 7. Non-perturbative improvement coeffi- 
cient CA- Results from one-loop bare and tadpole 
improved perturbation theory are denoted by dot- 
ted lines and crosses, respectively. 



which take into account the perturbative expan- 
sions MWi 



Zy = l- 0.1294305^ + 0(5^), 
Za = 1-O.116458go' + 0(5^). 



(52) 
(53) 



Eq. (|o|) and cq. ( [5l| ) are to be quoted with total 
errors of 0.5% and 1.0%, respectively. 

Fig.|^ shows the computed improvement coeffi- 
cient bv [pl and the fit [lT2|,f3l 



1 



r [y| and the fit 

- 0.7613. gg + 0.0012. 9^ - 0.1136 

1- 0.9145 5g ' ^ ' 



valid for < .go < 1- 

Unfortunately, 6a cannot be obtained in the 
same way, since the relevant Ward identity con- 
tains a physical mass dependence apart from the 
0(am) lattice artefact. We do not know how to 
separate the two. 

8. IMPROVEMENT OF THE VECTOR 
CURRENT 

For a complete determination of the renormal- 
ized improved vector current (p^, we further 
need the improvement coefficient cy. A compu- 
tation of cv is currently in progress p3|| . Here we 
outline the improvement condition employed and 
present preliminary results. 

The general idea is again quite simple. We have 
already seen above that correlation functions of 



0.1 


-0.1 
-0.2 ^ 
-0.3 
-0.4 



I I I I I I I I I I I 



I ' ' ' I 



preliminary 
I I I I I I I 







0.2 0.4 0.6 0.8 1 



Figure 8. Non-perturbative improvement coef- 
ficient Cv- One- loop bare perturbation theory is 
shown as a dotted line. 



axial current and vector current are related by 
chiral Ward identities. These are valid up to 
terms of order in the improved theory. Since 
the improved axial current is known, cy can be 
obtained through a suitable Ward identity. 

In order to excite states with vector quantum 
numbers, we use the boundary operator 



(0v)? = «^Ec(y)7feYC(z) 



(55) 



and define the bare, unimproved, correlation 
functions 



(56) 
(57) 



as well as the bare, improved correlation function 



fcAA(a;o,yo) 



18 ^ 

X 

x{iA,r,{x){A)i{y) {Oy)l) 



(58) 



These correlation functions are related through 
the Ward identity [ti < xq < 



(59) 

aa(^i>2;o)] + O(a^), 



Zy ky{xo) + acydQk^{xo) 
Za [fcAA(*2,a;o) - k\ 

which is valid in the massless limit. Eq. (^9|) is 
derived in much the same way as the Ward iden- 
tity discussed in the previous section. Inserting 
Za , Zy and ca , it can be solved for the improve- 
ment coefficient cy. 
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The reader will by now be aware that, for a 
precise definition of the improvement condition, 
we have to make specific choices for L,T,9, . . .. 
These, as well as a number of other details, will 
be given in ref. [Q. 

First results (cf. fig. ||) show cy to be much 
larger than the perturbative estimate eq. (|lj) for 
couphngs gQ > 0.9. Its effect on vector meson de- 
cay constants and semileptonic decay form factors 
may be significant in simulations with presently 
achievable resolutions. 

9. 0(a2) EFFECTS 

In sect. 2 we gave a perturbative example, 
where 0(a) improvement removed most of the 
cutoff effects, leaving only tiny effects of higher 
order in the lattice spacing. Of course, it is 
important to check in how far this is true non- 
perturbatively. 

9.1. The PC AC relation 

The first test of O(a^) effects is provided again 
by the PCAC relation. To set the scale, remem- 
ber that the cutoff effects in the PCAC mass m 
were as large as several tens of MeV before im- 
provement (cf. fig. The situation after im- 
provement, and for a somewhat larger value of 
the lattice spacing, is illustrated in fig. ^. Away 
from the boundaries, where the effect of states 
with energies of the order of the cutoff induces 
noticeable effects, m is independent of time to 
within ±2 MeV. 

A search for other lattice artefacts after im- 
provement yielded, as largest effect, the one de- 
scribed in the following. 

As chiral symmetry is violated by lattice arte- 
facts, there is no precise definition of a chiral 
point K = Kc = 1/(8-1- 2aTOc) for any finite value 
of the lattice spacing Rather, has an in- 
trinsic uncertainty, which is reduced from O(a^) 
to 0{a^) by non-perturbative improvement. 

Here, we define the value of k, where 

TO, eq. d^), vanishes tor C ^ C" ^ 9 = 0,T = 
2L,xo = T/2. To study the O(a^) effects we can 
then still vary the resolution a/L. 

Close to Kc the mass m is a linear function of 
K to a high degree of accuracy. The point Kc is 
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Figure 9. Unrenormalized current quark mass to 
in the improved theory, with non-perturbatively 
determined Csw and ca, as a function of the time 
xo on a 32 X 16^ lattice at /3 = 6.2 and k = 0.1350. 
The width of the corridor bounded by the dotted 
horizontal lines is 4 MeV. 



therefore easily found by linear interpolation (or, 
for /3 < 6.2, slight extrapolation). 

For (3 > 6.8, the results for a/L ~ 1/8 and 
a/L — 1/16 agree on the level of their statis- 
tical accuracy, which is better than 10~^. For 
lower (3, small but significant 0{a^) effects are 
seen (fig. They decrease rapidly as the lat- 
tice spacing is reduced.^ 

9.2. Hadronic observables 

It is of great interest to check the size of O(a^) 
terms in hadronic observables. At this conference 
first results for the improved theory have been 
presented by UKQCD |34| and we had the plea- 
sure to survey the calculations of Dc Divitiis et 
al. and Gockeler et al. |3^. Since all of these 
results are still preliminary, we restrict ourselves 
to some qualitative observations. 

Wherever a comparison is possible, the results 
of the different groups are in good agreement. As 
the calculations were initiated only fairly recently, 
results to date exist for two values of the lat- 
tice spacing only, which correspond to (a/rp)^ = 
0.035 and {a/rof = 0.019. Although the lattice 
artefacts should therefore vary by roughly a fac- 



^ Note that the dependence of 5kc on a/ro is non-trivial 
since there are two relevant nhysical scales, L and rg, in 
the regime covered by fig. |lOj. On general grounds one 



can only predict that Skc 
unknown function d. 



(a/ro) d{L/ro) with 
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0.5 



0.01 0.02 0.03 0.04 
(a/ro)= 

Figure 10. Change 5kc in the calculated value of 
Kc when L/a is increased from 8 to 16. The data 
have been obtained at /3 = 6.0, 6.2 and 6.4 (from 
right to left). For /3 > 6.8, 6kc is consistent with 
zero within the statistical errors that are at most 
10"^ 



tor of 2 between the two, a third, smaller value 
of a appears necessary to perform precise extrap- 
olations to the continuum limit. 

Nevertheless, one observes iQ,^ that the lat- 
tice spacing dependence of the ratio of the vector 
meson mass to the square root of the string ten- 
sion is reduced by improvement; indeed the val- 
ues calculated in the improved theory are close 
to those computed without improvement and ex- 
trapolated to a = 0. Unfortunately, it is not ob- 
vious that systematic errors in the determination 
of the string tension do not distort the picture. It 
would be desirable to repeat the analysis, replac- 
ing the string tension by p^ . 

The effect of the improvement of the axial cur- 
rent has been studied in the calculation of De- 
spite the rather small value of the coefficient ca, 
the admixture of represents a 10% change 
in /jr at (a/ro)'^ — 0.035. After inclusion of this 
effect, to/tt is independent of the lattice spacing 
within the statistical errors of around 3%. 

In summary no large lattice artefacts have so 
far been found in the (quenched) improved the- 
ory. Improvement (and of course proper renor- 
malization) of the composite fields is important 
in this respect. 



10. DYNAMICAL FERMIONS 

The ALPHA collaboration has started to im- 
plement improvement in full QCD with two 
flavours. 

In particular, we are currently calculating 
Csw(5o) along the lines described in section^ [Q. 
We simulate with the Hybrid Monte Carlo algo- 
rithm [ p7[ employing even-odd preconditioning of 
the Dirac operator and the Sexton- Weingarten 



scheme 38 1 for the integration of the equations 
of motion. The classical trajectories have unit 
length. For details of the implementation see 
ref . 1^ . Alternative algorithms are being studied 
at the same time jjo] ]. 

One slight difference from sect. || is that we do 
not insist on setting the improvement condition 
at M = 0, but rather keep aM small, say smaller 
than 0.02. Remember that the dependence of 
AM on M is insignificant (cf. fig. ||). The rea- 
son for not attempting to put M to zero exactly 
is not a difficulty in simulating in the massless 
limit. In the SF the massless Dirac operator has 
a gap [0, at least at not too large values of L, 
say L < 1 fm. Indeed, the simulations done so 
far show that it is irrelevant for the simulations 
whether we have aM = 0.02 or aM = 0.001. We 
do choose finite values of M in order to avoid an 
unnecessarily large effort for tuning M. 

Having said this, we mention, however, that as 
go is increased at fixed L/a = 8, the physical scale 
L increases and the gap in the Dirac operator be- 
comes smaller. This means that more Conjugate 
Gradient iterations are needed to obtain a solu- 
tion of the Dirac equation. Quantitatively, in one 
hour on a 256-node APE-100, we obtain around 
100 trajectories at = 0.8, 50 trajectories at 
go — I and "only" 30 trajectories at go = 1.1. 

At the same time our autocorrelation times for 
typical long-range observables grow from roughly 
1 trajectory at gQ = 0.5 to 4 trajectories at = 
1.05. 

The simulations were started at small values of 
go and reach gg = 1.1 by now. 

In 1-loop perturbation theory ||l^, there is a 
shift in the bare coupling /3 = 6 /go from the 
Wilson action to the improved action of —0.3 
{N{ = 2). Using this and estimates of the lattice 
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spacing with Wilson fermions |^], a rough esti- 
mate for the lattice spacing at = 1.1 is 0.1 fm. 
We emphasize that this is a rough estimate since 
the true renormalization may be quite different 
at such a large value of the coupling. Neverthe- 
less this indicates that Csw will not be needed for 
values of that are much larger than = 1.1. 

11. CONCLUSIONS 

We have been able to implement on-shcU 0(a) 
improvement non-perturbatively in (quenched) 
lattice QCD. The improvement coefficients Cgw, 
ca, the critical mass rric and the current nor- 
malization constants Zp^ and have been de- 
termined for bare couplings in the range < 
go < 1. In all cases, contact with bare pertur- 
bation theory could be made at couplings around 
(/q ~ 0.25 — 0.5. The convergence of perturbation 
theory appears to be significantly speeded up by 
tadpole improvement as shown by the crosses in 
figs. However, for values of /9 < 6.8, which is 
the range most relevant for present MC compu- 
tations, the quality of tadpole-improved pertur- 
bation theory is rather non-universal. There are 
examples such as 6v, ca, where perturbation the- 
ory is far off from our results. Renormalizations 
are large! 

This suggests that one should be cautious when 
estimating improvement terms to a low order 
of perturbation theory, as is currently done in 
several attempts to extend improvement to the 
0(a2) level. 

It is very pleasing that non-perturbative results 
can be obtained with statistical errors that are 
much smaller than the intrinsic O(a^) uncertain- 
ties of the perturbative results for the improve- 
ment coefficients. The elimination of these per- 
turbative uncertainties allows for reliable contin- 
uum extrapolations of the hadron spectrum and 
current matrix elements to be done in the future. 
Indeed, as a first step in this direction it was 
observed that residual O(a^) effects in hadronic 
quantities appear to be quite small. In particular, 
we point out that the a-independence of to/tt was 
observed after accounting for the relatively large 
correction proportional to ca, a term that would 
have been set to zero in tree-level improvement 



and is still tiny at 1 loop. 
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